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The energy loss of high pr partons propagating through a hot and dense medium is regarded as 
a valuable tool to probe the medium created in ultrarelativistic heavy-ion collisions. The angular 
correlation pattern of hadrons associated with a hard trigger in the region of pr ~ 1-2 GeV which 
exhibits a dip in the expected position of the away side jet has given rise to the idea that energy 
is lost predominantly to propagating collective modes ('Mach cones'). Recent measurements by the 
STAR collaboration have shown that for a high pt > 8 GeV trigger the angular pattern of associate 
hadrons for p-r > 4 GeV shows the emergence of the expected away side peak. These di-jet events 
suggest that the away side parton may emerge occasionally without substantial energy loss. Since 
in such a back-to-back configuration one of the partons may travel through the central region of the 
fireball, the average in-medium pathlength is substantial and the expected energy loss is not only 
sensitive to the initial geometry of matter but also to the change of geometry due to expansion. We 
show that radiative energy loss is able to explain the dijet events provided that the expansion of the 
medium is taken into account. 

PACS numbers: 25.75.-q 



I. INTRODUCTION 

Announcements have been made by all four detector col- 
laborations at RHIC 0] that a new state of matter, dis- 
tinct from ordinary hadronic matter has been created in 
ultrarelativistic heavy-ion collisions (URHIC). A new and 
exciting challenge for both experiment and theory is now 
to study its properties. The energy loss of hard partons 
created in the first moments of the collision has long been 
regarded a promising tool for this purpose 0,0,|3,|a,IE0- 
However, before this tool can be fully exploited, the in- 
fluence of the medium evolution on the energy loss has 
to be understood in a more quantitative way (cf. 
Recently, measurements of two-particle correlations in- 
volving one hard trigger particle have shown a surprising 
splitting of the away side peak for all centralities but 
peripheral collisions, qualitatively very different from a 
broadened away side peak observed in p-p or d-Au colli- 
sions . Interpretations in terms of energy lost to prop- 
agating colourless [13 in 111 and coloured J^J sound 
modes have been suggested for this phenomenon, and cal- 
culations within a dynamical model evolution have shown 
that the data can be reproduced under the assumption 
that a substantial amount of lost energy excites a sonic 
Shockwave [Tol| . 

However, the STAR collaboration has reported [l^ ITsI l 
that the picture of two-particle correlations is qualita- 
tively changed if both transverse momentum pT of the 
trigger hadron and transverse momentum cut of asso- 
ciate hadrons are increased. Specifically, for p^'^®""^ > 8 
GeV and p^'^"^ > 4 GeV a clear signal at the expected 
position of the away side jet has been observed. This sug- 
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gests that at such large momenta the away side parton 
is not always absorbed by the medium but may emerge 
with a substantial amount of energy left. Since radiative 
energy loss is based on a probabilistic picture of loosing 
the amount of energy AE' given the parton path through 
the medium, we note that there is always a finite (but 
small) probability that a parton does not loose any en- 
ergy even if it propagates through the densest central 
region. 

For a realistic theoretical model of dijet events it is nec- 
essary to note that for a typical back-to-back configu- 
ration the length of propagation through the medium 
can be substantial: If the hard vertex is close to the 
surface and the near side parton travels outward (as is 
the case for most triggered hadrons), the away side par- 
ton has to cross the whole expanse of the medium, i.e. 
a length of ^ 10 fm before it reaches the opposite sur- 
face. At such length- and timescales, the problem cannot 
be approximated by a static situation: both longitudinal 
and transverse flow substantially alter the distribution of 
matter during the propagation. To illustrate the point 
more clearly: Naively one would expect that such long 
pathlengths lead to extremely high energy loss due to 
the quadratic pathlength dependence of the energy loss 
on the pathlength L. However, due to the accelerated 
longitudinal and transverse expansion, the medium en- 
ergy density drops with a high power of time (and hence 
L), more than cancelling the quadratic pathlength de- 
pendence. Thus, in a dynamic picture penetration of the 
fireball core region is possible though rare. 

In the following, we present the formalism used to calcu- 
late the expected high px angular correlation pattern in 
the presence of an evolving medium and compare with 
the experimental results. 
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II. PARTONIC ENERGY LOSS 

We aim to describe data involving trigger hadrons beyond 
Pt = 8 GeV and associate hadrons with > 4 GeV. 
We assume that at such energies the observed 2-particle 
correlations are dominated by three contributions: 1) an 
away side parton emerging from the medium with a large 
fraction of its original energy left undergoing fragmen- 
tation ('punchthrough') 2) a subleading hadron created 
in the fragmentation of the trigger hard parton ('near 
side associate') and 3) a subleading hadron created in 
the fragmentation of the away side parton ('away side 
associate', only for punchthrough).). We posit that the 
effect of the medium can completely be treated on the 
partonic level by inducing partonic energy loss but does 
not substantially influence the fragmentation process. 
We note that any parton fulfilling the trigger and lead- 
ing to associated yield must have an energy in excess 
of 12 GeV (typically 15-20 GeV). Making a simple es- 
timate for the length scale for fragmentation as Lfrag = 
c • Epart I Qhad ^i^'^ Qhad ~ 1 GcV wc find length scales of 
~ 4 fm which is enough to carry most partons out of the 
dense regions of the fireball, thus the assumption seems 
not unreasonable. 

On the other hand, based on the low correlation pat- 
tern and our results in ^3 we argue that the dominant 
part of the energy lost to the medium does not appear as 
an accompanying cone but rather excites collective modes 
of the medium and hence is observed in the momentum 
region below 1.5 — 2 GeV where the medium shows still 
collevtivity (indicated by the validity of a hydrodynam- 
ical description for the spectra, cf. e.g. [13). Thus, 
we do not expect the 'lost' energy to contribute substan- 
tially to high Pt yield in addition to the sources 1) to 3) 
listed above. 

For the time being we focus on central collisions only. 
Key quantity for the calculation of jet energy loss is 
the local transport coefficient q(r]s,r,T) which charac- 
terizes the squared average momentum transfer from the 
medium to the hard parton per unit pathlcngth. Since wc 
consider a time-dependent inhomogeneous medium, this 
quantity depends on spacetime rapidity rjs = ^ In , 

radius r and proper time r = \Jt'^ — (for the time be- 
ing we focus on central collisions only). The transport 



coefficient is related to the energy density of the medium 
as g = ce'^/'* with c = 2 for an ideal quark-gluon plasma 
(QGP) 01 . In the present study, we have adjusted c such 
that the observed nuclear suppression is reproduced, re- 
sulting in c = 4. 

In order to find the probability for a hard parton P(A£') 
to lose the energy Ai? while traversing the medium, we 
make use of a scaling law ^ISj which allows to relate the 
dynamical scenario a static equivalent one by calculating 
the following quantities averaged over the jet trajectory 
C(r): 
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as a function of the jet production vertex Tq and its angu- 
lar orientation (p. We set q = whenever the decoupling 
temperature of the medium T = Tp is reached. In the 
presence of flow, we follow the prescription outlined in 
jsLITgll and replace 

g = ce3/4(p)^ce3/4(T"^»^) (3) 

with 

=p(e) + [e+p(e)]^^ (4) 

where (3± is the spatial component of the flow field or- 
thogonal to the parton trajectory. In the above two ex- 
pressions, the spacetime dependence (77^ , r, r) of pressure 
p and energy-momentum tensor r"^"^ have been sup- 
pressed for clarity. 

Using the results of |23|, we obtain P{AE) from ujc and 
(qL) as a function of jet production vertex and the angle 
(j) from 
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which makes use of the distribution Loj^oi gluons emit- 
ted into the jet cone. The explicit expression of this 
quantity for the case of multiple soft scattering can be 
found in 



The medium enters via the local energy density e(77s, r, t) 
in Eq. and flow j3±_ transverse to the jet propagation 
direction in Eq. and it is immediately obvious from 
Eqs. imni that a strongly dropping energy density as a 
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function of space or time will cancel the quadratic path- 
length dependence of e.g. lo^ in the case of a constant 
9- 

To describe the medium we use a parametrization of the 
spacetime evolution which is based on a simultaneous fit 
to hadronic spectra and Hanbury-Brown-Twiss correla- 
tion radii at low '2l'| (see also for details of the 
formalism) and which has been used to successfully pre- 
dict the emission of thermal photons [2^ . We have 
used the same description to study the excitation of Mach 
cones by high pT jets IQj and the sensitivity of the nu- 
clear modification factor 



momenta of the colliding partons i and j are given by 



RAA(pT,y) 



d^N^^/dprdy 
TAAiO)d^(r^^/dpTdy 



(6) 



to assumptions about flow in the evolution dynamics Q ■ 
We therefore believe that the model gives a fair represen- 
tation of the relevant physics of bulk matter in both the 
early and late stages of the evolution and incorporates 
the correct amount of quenching demanded by Raa- 



III. MONTE CARLO SAMPLING OF THE 
TRIGGER CONDITIONS 



The approach chosen here is similar to [23 in the sense 
that we perform a Monte Carlo (MC) sampling of the 
trigger conditions; with the crucial difference that the 
model presented here includes the full transverse and lon- 
gitudinal expansion dynamics. 

For central collisions, we assume without loss of general- 
ity that the near side parton always propagates into the 
negative x direction in the transverse {x,y) plane. If we 
define the —x direction as tt, the away side parton then 
propagates almost at angle zero, modulo the intrinsic fcy 
of order ~ 1 GeV (the resulting angular spread is small 
given that we're mainly interested in partons of 15 GeV 
and above). 

We start by creating jet vertices based on the nuclear 
overlap function TaaO^) — J dzp'^{h, z). Next we spec- 
ify the jet energy Ej^t by randomly sampling the pQCD 
parton spectra above the trigger energy. 
Basis for the sampling is the leading order (LO) pQCD 
expression for the back-to-back production of two partons 
in the transverse plane 



da 



AB^kl+X 



dp^dyidy2 



^Xif,/A{xi)x2fj/B{x2)^{s,i,u). 



dt 



(7) 

where fi(j)/A(B)i^i(j)j Q^) denote the nuclear parton dis- 
tributions 28] (in which we have suppressed the scale de- 
pendence on for clarity in the above expression) and 

(s, i, it) are the leading order pQCD subprocesses 
for all different combination of parton types i, j, fc, I (ex- 
plicit expressions can be found in [2^). The fractional 
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xi,2 = —r= [exp(±2/i) 
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In "271 it has been shown that this framework, when 
folded with the fragmentation functions and supple- 
mented with a if— factor to adjust overall normaliza- 
tion, is able to describe inclusive production of charged 
hadrons in p-p and p-A collisions for a wide range of en- 
ergies. This indicates that NLO effects can largely be 
absorbed into the if —factor. Note that such a constant 
normalization factor drops out in ratios such as Raa or 
the yield per trigger of dijets we're interested in, thus LO 
pQCD can be expected to work reasonably well. How- 
ever, we do introduce intrinsic kx to account for the cor- 
rect geometry and the broadening of the observed away 
side jet (see below). 

From Eq. (|7J) we also determine the probability that the 
parton is a gluon (by setting e.g. k = g). We propagate 
the near side parton to the surface and determine uJc and 
{qL) by evaluating Eqs. along the path. The result- 

ing values serve as input for the probability distribution 
of energy loss P{AE) as determined in j23|. Note that 
in addition to a continuous distribution P(AE) there is 
also a discrete part P{Q) which reflects the probability to 
loose no energy. This part is finite (but often small) for 
all possible vertices and paths. 

Often the plasma frequency lUc is far above the avail- 
able jet energy and P{AE) thus extends to energies 
AE ^ Ejct- This reflexts the fact that the radiative en- 
ergy loss in ^0 ] is derived in the limit of infinite parton 
energy and translates into an uncertainty. Two differ- 
ent prescriptions have been considered as the upper and 
lower limit quantifying this uncertainty: 

• Reweighting: Here, one truncates P{AE) at AE — 
Ejct and renormalize to unity by dividing out the 
factor P{AE). 

• No-reweighting: Alternatively, one may truncate 
P{AE) at AE = Ejct and consider the parton en- 
ergy to be zero if AE exceeds Ejct ■ 

The quenching power of the medium is higher in the non- 
reweighted case. However, note that reweighting is only a 
reasonable prescription if Ej^t is above the AE at which 
P{AE) peaks. If that is not the case (as regularly ob- 
served for away side partons propagating through dense 
regions in the simulation), reweighting becomes mean- 
ingless as easily seen from the fact that under some con- 
ditions the average energy loss after reweighting may de- 
crease if q is increased. On the other hand, we are not 
interested in details of energy losses of order 0{Ejet) - 
once the energy loss reaches a sizeable fraction of the par- 
ton energy, the event will fall below the trigger/associate 
threshold and the actual AE is irrelevant. Note that this 
is very different from a calculation of Raa where the par- 
tons experiencing substantial energy loss do not fall be- 
low a threshold but reappear elsewhere in the spectrum. 
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The dijet measurement is only sensitive to the probabil- 
ity to have little or no energy loss, and this should be 
independent of reweighting considerations as AE is al- 
ways a small fraction of Ej^t for almost all hadrons above 
the trigger or associate cut. Thus, we do not apply the 
reweighting in the following. 

We determine the actual energy loss of the near side par- 
ton by sampling P{AE). To find the energy of the lead- 
ing hadron, we use the KKP fragmentation functions 
to find the probability Di{z, /i) for a parton i to fragment 
into a hadron carrying the momentum fraction z at a typ- 
ical hadronic scale fi. Thus, we evaluate the schematical 
expression 

/ 

(9) 

in a MC framework where da^^~*^^^ is the pQCD ex- 
pression for the hard scattering process. 
If the resulting hadronic Phad = ^Pjot ~ zEj^t fulfills the 
trigger condition we accept the event and proceed with 
the calculation of associated hadrons and the away side 
parton, otherwise we reject the event and continue the 
MC sampling by generating a new vertex. 

IV. MONTE CARLO SAMPLING OF THE 
AWAY SIDE AND ASSOCIATED HADRONS 

If an event fulfiUig the trigger has been created, we deter- 
mine the intrinsic kx being added to the away side parton 
momentum. We sample a Gaussian distribution chosen 
such that the widening of the away side cone without a 
medium is reproduced. Since this is a number of order 
1 GeV whereas partons fulfilling trigger conditions have 
frequently in excess of 15 GeV we note that this is a small 
correction. 

We treat the far side parton exactly like the near side par- 
ton, i.e. we evaluate Eqs. |2J) along the path and find 
the actual energy loss from P{AE) with ujc, (qL) as in- 
put. If the away side parton emerges with a finite energy, 
we use the fragmentation function Di{z, /i) to determine 
the momentum of the away side hadron. If this momen- 
tum fulfills the cut for associated particle production, we 
count the event as 'punchthrough'. 

In addition, we allow for the possibility that the frag- 
mentation of near and away side parton produces more 
than one hard hadron. Since we are predominantly in- 
terested in the quenching properties of the medium, we 
sample this probability using the measured probability 
distribution Ajizp) of associated hadron production in 
d-Au collisions [3 ^| Eis a function of zt where zt 
is the fraction of the trigger hadron momentum car- 
ried by the associated hadron. We include a factor 
0{Ejct — -Etriggcr — AE — Sassoc) On the near side and 
d{Ejct — -E-punch — AE — -Eassoc) On the far side to make 
sure that energy is conserved. Note that associated pro- 
duction on the far side above the pT cut is only possible 



if a punchthrough occurs. We count these events as 'near 
side associate production' and 'away side associated pro- 
duction'. 

Both the fragmentation function Di{z,fi) and the as- 
sociated production probability distribution Ai{zp) are 
steeply falling functions of z, zp. They have to be regu- 
lated by a lower cutoff, otherwise the fragmentation func- 
tion will be evaluated in a regime where it is not applica- 
ble (In calculations starting from the hadron spectrum, 
a lower cutoff at the hadronic scale Pt, i.e. Zmin{PT) = 
arises from the fact that no parton can exceed the pr 
of the kinematic limit. In the present Monte Carlo model 
where the parton spectrum is determined first it is diffi- 
cult to implement the same hadronic Pt dependence of 
the cutoff, hence we follow a more heuristic approach.). 
We choose the cutoff in such a way that the measured 
yield per trigger as a function of the cut on associate 
particle production is reproduced in the case of d-Au col- 
lisions, i.e. without energy loss. We find that z""" = 0.05 
and z™'" = 0.1 gives a fair representation of the data and 
use the same values also in the presence of a medium. 
The resulting difference then only reflects the effect of 
the medium on the parton energies. 
Thus, the yield per trigger on the near side is determined 
by the sum of all 'near side associate production' events 
divided by the number of events fulfilling the trigger, the 
yield per trigger on the away side is given by the sum 
of 'punchthrough' and 'away side associated production' 
divided by the number of events. These quantities can 
be directly compared to experiment. 

V. RESULTS 

In Fig. n we show the resulting distribution of vertices 
leading to a valid event and the position of vertices for 
which a dijet was observed. We can infer a few interesting 
observation from the distribution: 

First, most triggered events cluster around the near side 
surface within a 'skin' of about 3 fm thickness. This is 
expected, as this is the region in which the energy loss of 
the near side parton is expected to be small. However, a 
sizeable number of triggered events also originates from 
the deep central region. This 'skin' effect has already 
been discussed in e.g. psj . 

More interesting is the distribution of events leading to 
an associated away side hadron ('punchthrough'). This 
distribution is manifestly different from the vertex dis- 
tribution for triggered events. First, its main strength 
comes from the region around a; ~ 0, hence events with 
equal pathlengths on the near and away side are more 
likely to result in dijets. Second, it extends much further 
into the large R region (the 'halo') than the distribution 
of triggered vertices. In the halo (i.e. in the tail of the 
nuclear overlap distribution the energy density is low to 
start with, and by the time a parton from the halo prop- 
agates into the center the energy density there has been 
reduced sizeably by the fireball expansion, hence such 
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Triggered vertex distribution Punchtlirougli distribution 
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FIG. 1: Left panel: Normalized density of vertices for events with a trigger hadron above 8 GeV. Right panel: Normalized 
density of vertices for which an associate hadron in the momentum region 4 < pr < GeV is found. The densities are based on 
8600 simulated events, in all cases the near side parton propagates into negative x direction. 



a parton never encounteres significant energy densities 
and dijet events are hence favoured in this region. In 
addition, there may be some tendency for dijet events to 
accumulate in the surface region close to a; = 0, y ±6 
fm where the pathlength through matter is small for both 
near and away side parton. Note that while both plots 
are in principle expected to be symmetrical around the 
X-axis, in practice they are not due to the finite size of 
the sample of simulated events. 

Let us briefly remark on the mechanism which under- 
lies this distribution: We observe that about 75% of all 
trigger partons probe the discrete part of the quenching 
weight, i.e. did not undergo any energy loss. This is in 
agreement with the observation made in psj that the en- 
ergy loss probability of observed hadrons is dominated by 
the discrete part. However, the crucial difference is that 
the expansion dynamics and subsequent drop in density 
allows the region in which the discrete quenching weight 
is still relevant to move further into the fireball center. 

This in turn has a strong influence on the pathlength dis- 
tribution of the away side parton — if the production ver- 
tex is not confined to the surface, the average pathlength 
of the away side parton is significantly reduced (and in- 
deed the distribution favours about equal pathlengths 
for near and away side parton). The magnitude of en- 
ergy loss on the away side for observed hadrons is rather 
strongly momentum dependent: While only ~ 30% of 
all hadrons in the 4-6 GeV momentum bin probe the 
discrete quenching weight in the simulation, the fraction 
increases to ^ 60% in the 6-1- GeV associate hadron mo- 
mentum region. Thus, the larger the difference between 
trigger energy and associate hadron cut, the more is the 
sensitivity to the continuous part of the energy loss prob- 



ability increased. 

We show the resulting yield per trigger compared with 
the STAR data in Fig. |21 (since the model calculation 
has only been performed for central collisions and van- 
ishing medium, we cannot address the full range of Npart 
probed by the experiment). The results of the model 
calculation agree surprisingly well with the data, consid- 
ering that all hadronization parameters have been fixed 
without reference to the medium and that all quenching 
properties of the medium are completely determined by 
the description of Raa 

It is reassuring that the model works especially well for 
the pt > & GeV associate hadron binning, as in this 
region the assumption of fragmentation being the only 
relevant mechanism for hadronization is better justified. 
It is plausible that the discrepancy of ~ 30% in the 4-6 
GeV window on the away side can be traced back to a 
recombination contribution (not contained in the present 
model) as the hadron-species averaged contribution of 
recombination to the hadron spectrum in central Au-Au 
collisions bewteen 4 and 6 GeV is of the same order of 
magnitude (cL Fig. 4 in 3_Q] ) . 



VI. CONCLUSIONS 

We have evaluated radiative energy loss for high pT par- 
tons within a dynamical model for the evolution of hot 
and dense matter. The same energy loss formalism which 
describes Raa is in such a model also able to describe the 
correct amount of experimentally observed high px trig- 
gered high pt associated hadron ('dijet') events under 
the assumption that these events represent hard partons 
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FIG. 2: Yield per trigger on the near side (left panel) and away side (right panel) for trigger hadron above 8 GeV in the model 
calculation as compared with STAR data. 



emerging from the medium and fragmenting in vacuum. 
Note that in this picture, the central region is still on 
average very opaque. The average value of the energy 
loss for a quark propagating outward from the center of 
the fireball is about (Ai?) = 23 GeV - on average partons 
from the center cannot reach the surface. Thus, most 
of the observed trigger and dihadron correlation yield 
coming from the central region of the medium represents 
a relatively rare class of events in which the low Ai? tail 
of the probability distribution or the discrete quenching 
weight is probed. 

We have argued based on the distribution of produc- 
tion vertices that a sizeable amount of observed asso- 
ciated hard hadrons propagates a sizeable distance > 5 
GeV through the medium. For such long timescales, the 
drop in energy density from both longitudinal and trans- 
verse expansion cannot be neglected, therefore a dynam- 
ical model evolution is crucial. Especially events leading 
from the halo region which never encounter a significant 
energy density do not emerge in a static description of 
jet quenching. Another crucial ingredient of the simu- 
lation is the discrete part of the quenching weight, i.e. 
the probability to have no energy loss in spite of a size- 
able pathlcngth through the medium. The magnitude of 
this probability likewise depends on the line-integral over 
the dropping medium density and is thus also on aver- 
age significantly larger in a dynamical picture than in a 
static model. While this probability is still numerically 
small, in the calculation of dijet yields it competes with 
other small quantities which allow for hard dihadron pro- 
duction like the probability to select a parton from the 
high j)T tail of the parton spectrum or to fragment into a 
hadron with z « 1, probing the tail of the fragmentation 
function or to select a vertex at large radius where Taa is 
small. Since the dihadron yield itself is not numerically 
large, all these different effects contribute in a way that 
is characteristic for the medium evolution. 
We conclude that radiative energy loss is well able to 
account for both Kaa and the observed pattern of dijet 



events. While the central part of the medium is very 
dense, the model simulation indicates that it is by no 
means completely black, and high triggered events 
seem indeed to probe the densest part of the medium to 
some degree. 

Recently, elastic energy loss has been suggested as an 
additional ingredient to describe the nuclear suppression 
(see e.g. [sj)- To the degree that the magnitude of this 
contribution is parametrically given by AE' ~ d£,q{£,) 
(note the absence of a factor f as compared with Eq. 
for radiative energy loss), the results of the present inves- 
tigation place tight limits on the importance of coUisional 
energy loss: Since the yield per trigger of hard dihadrons 
is determined by a delicate interplay between differences 
in average near and away side pathlength, in turn leading 
to different energy loss probabilities, the scenario cannot 
easily accomodate a different pathlength dependence. In 
particular, the coUisional integral above is dominated by 
early times when the density is largest (whereas the ra- 
diative integral Eq. ^ gets most support around 3-4 
fm/c) and is completely insensitive to late times (and 
consequently long paths). However, this implies that 
there would be no additional suppression for the away 
side partons in dihadron measurements as compared to 
d-Au; while the total yield would go down like Raa, the 
yield per trigger (which in the present simulation rests 
on the systematic difference in pathlength for near and 
away side parton and its weighting via Eq. (^) would 
be unchanged. A more detailed future investigation will 
make this point more quantitative. 
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